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Abstract 

We compute the 0(a s ) QCD corrections to the partonic process gg — * tiZ at the LHC. This 
partonic channel is the dominant component of the scattering process pp — > ttZ, which will be 
important for measuring the ttZ electroweak couplings. The 0{a s ) corrections increase the total 
cross section by up to 75% for reasonable choices of the renormalization and factorization scales. 
Inclusion of these contributions descreases the residual scale dependence of the cross section coming 
from uncalculated higher order terms to ±5%. 
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I. INTRODUCTION 



The discovery of the top quark in experiments at the Fermilab Tevatron established 
the validity of the Standard Model as a quantum field theory based on the gauge group 
SU(3) x SU(2) x U(l) with three generations of quarks [lj]. Subsequently, the production 
cross section for ti pairs in proton anti-proton collisions was measured with good precision 
and the value of the top quark mass was established to nearly the percent level [2J. Very 
recently, the first observation of single top production was reported [3|. 

While the above studies are important milestones in top quark physics, these measure- 
ments probe only a limited subset of top quark properties. The process pp — > ti at the Teva- 
tron is sensitive to the top quark mass and color charge. The decay t — > Wb is marginally 
sensitive to the tWb coupling. Other important properties of the top quark, such as its cou- 
plings to the Z-boson and the photon and its Yukawa coupling, cannot currently be directly 
accessed experimentally. These couplings are sensitive probes of new physics effects; in par- 
ticular, the Zti coupling is affected by tree-level mixing with additional Z gauge bosons 
and vector-like fermions. 

While the best place to perform detailed studies of the top quark is the International 
Linear Collider, useful information about the top quark couplings to electroweak gauge 
bosons and the Higgs boson can be obtained from the LHC where the processes pp — > tiH, 
pp — > tiZ and pp — > tij can be observed. However, in order to use those observations, 
accurate theoretical predictions for the signal and background events are required. Since 
the cross sections for all signal processes scale as cr LO ~ a 2 s at leading order in the QCD 
perturbative expansion, there is significant uncertainty in the normalization of the leading 
order cross section. To remove this theoretical error, a computation of the cross section 
through next-to-leading order (NLO) in a s is required. 

Recently, the possibility of measuring the top quark couplings to the photon and the 
Z-boson at the LHC was studied in detail j3, 0] . The study in Ref . 0] required that the Z- 
boson decays leptonically, and allowed for either zero or one leptonic W boson decay. Both 
channels have a reasonable signal to background ratio. It was pointed out in Ref. [5( that 
further improvements can be achieved by considering the decay Z — > vv. The background 
processes to pp — > tiZ and pp — > tij scale with a high power of a s , Obckg ~ a™, n > 2. 
Consequently, theoreticalpredictions for these processes are quite uncertain. Nevertheless, 
as pointed out in Refs. |4|, |5|], it may be possible to experimentally study them in signal- 
free regions to control the background normalization. In principle, this technique allows 
the extraction of the tiZ and tij couplings from the tiZi^j) cross section measurements, if 
not for the normalization uncertainty in the signal cross sections themselves associated with 
unknown QCD corrections. It was shown in Ref. 0] that this uncertainty is a limiting factor 
in extracting anomalous tiZ couplings in the leptonic Z decay channels. 

In this paper we take the first step towards computing the NLO QCD corrections to tiZ 
production at the LHC by considering the gg — > tiZ partonic subprocess. We focus on 
this reaction for two reasons. First, gluon collisions at the LHC are important because of 
the large gluon luminosity. This can be seen already at leading order where the gg — > tiZ 
subprocess gives approximately 60% of the full result. Second, the computation of NLO QCD 
corrections to the gg — > tiZ subprocess is technically more involved than the corresponding 
computation for quark anti-quark annihilation channels. It can also be expected that NLO 
QCD corrections to gg — > tiZ are larger than corrections to qq — > tiZ, since the color charge 
of gluons is larger than that of quarks. For these reasons we believe that the computation 
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presented in this paper gives a good idea about the importance of 0(a s ) corrections to the 
full pp — > tiZ process. However, a detailed study of how NLO QCD influences the extraction 
of the tiZ coupling requires the inclusion of the qq and qg channels. We plan to extend our 
computation to incorporate the quark-initiated production processes in the near future. 

Our paper is organized as follows. In the next Section we introduce our notation and 
briefly discuss the calculation of the leading order cross section. Section III describes the 
computation of the 0{oc s ) corrections. In Section IV, numerical results are presented and 
discussed. We conclude in Section V. 



II. NOTATION AND THE LEADING ORDER CROSS SECTION 

We consider the process p(P\) + p{P2) ~^ ttZ . The factorization theorems for hard 
scattering processes in QCD allow us to write 



da = / dx 1 dx2fi{x 1 )ff{x 2 )da ij _> t i Z {x 1 X2S), 

ij o 



(1) 



where f-J (xi j2 ) are the parton densities that give the probability to find parton in 
the proton 1(2) with momentum p^ = £i(2)-Pl(2)- The center-of-mass energy squared of the 
proton-proton collision is introduced in Eq. (TjQ), S = 1P\ ■ P 2 . In this paper, we only consider 
gluon collisions at the LHC, so we set i — j — gluon in Eq. ([T]) and use VS = 14 TeV. 
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FIG. 1: Sample diagrams for gg — > ttZ appearing at leading order. 



At leading order in the a s expansion, there are eight Feynman diagrams that contribute 
to gg — > tiZ; some examples are shown in Fig.1. The computation of the leading order cross 
section is straightforward. We use QGRAF [6( to generate the relevant Feynman diagrams 
and then MAPLE and FORM [7] to manipulate this output. Throughout the paper we set 
m t = 170.9 GeV, m z = 91.18 GeV, and raw = 80.45 GeV. For the coupling of the Z-boson 
to quarks, we employ 



Zqq 



lm 2 w G F 



v^cos 2 ^ 



{9l + 9 q a J 5 ) 



(2) 



where gl 



rpq 

J 3 



— Q q sin 2 9w, g q a = — sin 2 9yy = 1 — m^/ml = 0.2215 is the sine squared 



of the electroweak mixing angle, Tg is the weak isospin of the quark q, Q q is the electric 
charge of the quark q in units of the proton charge and Gp is the Fermi constant. Numerical 
results for the leading order cross section are reported in Section IV. 
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III. NEXT-TO-LEADING ORDER COMPUTATION 

We now discuss the computation of the next-to-leading order corrections to the partonic 
process gg — > tiZ. Three distinct contributions should be considered: virtual corrections to 
the leading order process gg — ► ttZ, the real emission corrections gg — ► tiZ + g, and the 
renormalization of the leading order cross section. 

We first consider the computation of the virtual corrections to gg — > tiZ. There are 
162 diagrams that contribute at next-to-leading order. Some examples of these one-loop 
diagrams are shown in Fig. 2. These diagrams fall into two categories: 

• diagrams in which the Z-boson couples to the top quark, such as the first two shown 
in Fig. 2; 

• diagrams in which the Z-boson couples to an internal quark propagator, such as the 
final one in Fig. 2. 

The first class contains contributions proportional to the coupling constant combinations 
{gl) 2 and {gi) 2 - Diagrams in this class can be computed with standard (^-dimensional Dirac 
algebra, using the anti-commutation prescription for 75. The second class contains contribu- 
tions proportional to g\g% and g\g% where q denotes any of the six quarks. The diagrams in 
this class contain products of two Dirac traces with single 75 terms, leading to contributions 
containing a product of two Levi-Civita tensors. Although this second set of diagrams is 
finite, we use the <i-dimensional prescription for 75 described in Ref. [§[ for calculational 
convenience in the intermediate steps. 




FIG. 2: Sample diagrams for the 0{a s ) virtual corrections to gg — > ttZ. 

To compute the virtual corrections, we employ the method described in Ref. [9|. We sum- 
marize here the salient features of this technique. Further details can be found in Ref. [||. 
We compute the square of the scattering amplitude, summing over the fermion spins and 
polarizations of gluons and the Z-boson. For each one-loop diagram interfered with the 
full Born amplitude, we perform the integration over the loop momentum and arrive at 
an integral over Feynman parameters. The integration over Feynman parameters cannot 
be computed numerically because it diverges. While the ultraviolet divergencies factorize 
after the integration over the loop momentum and therefore have no effect on the integra- 
tion over Feynman parameters, the infrared and collinear singularities appear when certain 
Feynman parameters approach boundary values. To have a representation of the integrals 
suitable for numerical evaluation, we must extract the infrared and collinear singularities. 
We accomplish this using the method of sector decomposition |Io| . However, even after 
these divergencies are extracted, it is not possible to perform the Feynman parameter in- 
tegrations numerically because there are singularities inside the integration region coming 
from internal loop thresholds. These singularities are avoided by deforming the integration 
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contour into the complex plane following the suggestion in Ref. [11]. Once infrared and 
collinear singularities are extracted and the integration contour is deformed, we obtain rep- 
resentations of the Feynman parametric integrals suitable for numerical integration. This 
integration is performed using the adaptive Monte Carlo integration algorithm VEGAS [l2| 
as implemented in the CUBA library [131 ]. 

We briefly comment here on the stability of the numerical integration. In principle, after 
extracting infra-red and collinear singularities and deforming the contour into the complex 
plane, the integration should be numerically stable. For most phase space points that we 
consider, we find this to be the case. However, there are situations where straightforward 
application of the algorithm described in 0] leads to numerical instabilities. These instabil- 
ities typically occur when the sector decomposition procedure "overdoes" the extraction of 
singularities, splitting the integral into too many sectors which each contain unphysical sin- 
gularities. Problems of these type occur typically in the evaluation of five-point functions. 
We handle these problematic cases by choosing a more sophisticated parameterization of 
the original Feynman integral. For example, if in a five-point function two propagators Di 2 
never become singular, it is beneficial to combine them first into a single propagator: 

l 

<U D® = + -£)D 2 . (3) 



D l D 2 J D(0 2 
o 

Using this and similar techniques we are able to obtain a stable numerical representation 
for all virtual diagrams that contribute to the gg — > tiZ process. 

In order to arrive at a finite result at next-to-leading order, renormalization is required. 
We employ the MS renormalization constants for the QCD coupling a s with five active 
flavors; the top quark contribution to the coupling constant renormalization is decoupled 
through a zero-momentum subtraction. In addition, we renormalize the top quark wave 
function and the top quark mass on shell. No other renormalization constants or counter- 
terms are required to render 0(a s ) corrections to gg — * ttZ process ultraviolet finite. 

We must also compute the real emission diagrams corresponding to the process gg — > 
tiZ + g. There are 50 such diagrams. When the gluon in the final state is soft or collinear to 
the gluons in the initial state, the matrix element for gg — > ttZ+g process becomes singular. 
To compute the corresponding contribution to the cross section, we split the phase-space for 
real gluon emission into soft, collinear and hard parts, using the two cutoff slicing method 



14J . The singularities that originate from collinear gluon emission are removed using the 
MS renormalization of the parton distribution functions. The computation of real emission 
corrections to gg — > ttZ using the two cutoff slicing method is analogous to the computation 
of real emission corrections for the gg — > tiH process described in [15J, where many details 
can be found. 



IV. RESULTS 

We now present the results of our computation. For all numerical results reported in 
this paper, we set the renormalization scale \i T and the factorization scale fif equal to a 
common scale: \i r = /// = //. We employ the MRST parton distribution functions [l6[ at 
the appropriate order in the perturbative expansion. We have applied a number of checks 
to our calculation. 
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1. We have compared the leading order cross-section obtained with our code with the 
result of a similar computation using the program MadEvent [17j] and have found 
complete agreement. 

2. The divergent parts of the NLO virtual correction to gg — > ttZ can be expressed in a 
simple analytic form, as explained in detail in We have checked that our results 
satisfy this check. 

3. We have checked that all divergences cancel once the real emission processes, the 
collinear counterterms, the renormalization constants, and the virtual corrections are 
combined. 

4. An important check of our NLO virtual result is provided by its independence of the 
size of the contour deformation. 

5. Finally, we have implemented all parts of the computation in several different codes 
that agree for all observables studied. 
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TABLE I: The cross section at leading and next-to-leading order for various values of the renor- 
malization and factorization scales. The numerical error on all numbers presented is 1% or better. 



We begin with the total cross section. The results of our calculation are summarized in 
Table [H We show there the total cross section computed through leading and next-to- leading 
order for the five values of the scale fj, — (1/8, 1/4, 1/2, 1, 2)fi , where /i = (2m t + m z ). For 
the QCD coupling constant, we use values consistent with the MRST parton densities at 



leading and next-to-leading order |16fl at fi = mz and evolve these values using the beta- 
function at the appropriate order in the perturbative expansion. We see from Table Hthat the 
leading order result changes by more than a factor of 3 when the scale \x changes between its 
minimal and maximal values. This suggests significant 0{pt s ) corrections. Indeed, as clearly 
seen from the last column in Tabled the NLO QCD corrections are large. They change the 
leading order corrections by up to a hundred percent, depending on the value of the scale 
fi. The dependence of the next-to-leading order cross section on the scale \x is significantly 
reduced; when \i is varied between fio/8 and 2/z , cr NL o changes by only a factor of 1.25. 

While formally all values of the renormalization and factorization scales are allowed, the 
best predictions for physical quantities at low orders in perturbation theory are obtained 
when the scales are set to values suggested by the physics of the problem. The typical 
partonic center-of-mass energies leading to the tiZ final state are moderate, ^/s^ ^ 10 3 GeV, 
and the typical transverse momenta for the Z-boson and top quarks are 100 — 200 GeV. 
It seems reasonable to choose scales comparable to the typical transverse momenta. This 
suggests setting \i ~ few x 100 GeV. Since [Aq ~ 2m t + mz ~ 450 GeV, we consider 
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[i = [/xo/4, //o] as a sensible range for the variation of the renormalization and factorization 
scales. As follows from Table D], in this range o"lo(cnlo) changes by about ±30%(±5%) if 
measured from the central value /i = /xq/2. 
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FIG. 3: Left panel: the transverse momentum distribution of the Z-boson at leading or- 
der and next-to- leading order in a s . Right panel: the dependence of the K- factor defined as 
K = do"NLo/dc7Lo on the transverse momentum of the Z-boson. We set fj, = fJ,o/2 in both plots. 

As pointed out in Ref . 0, 0] , the signal-to-background ratio for pp — > tiZ production is 
improved if a cut on the transverse momentum of the Z-boson p±(Z) is imposed since many 
background processes decrease faster than the signal when p±{Z) increases. It is therefore 
interesting to study the effect of the NLO QCD corrections on the p±(Z) distribution. In 
the left panel of Fig. 3 we show the transverse momentum distribution of the Z-boson 
at leading order and next-to- leading order in the strong coupling constant. We combine 
events in 20 GeV bins. We set the factorization and renormalization scales to /x = /j.q/2. 
While the QCD corrections change the distribution significantly, they affect only the overall 
normalization. This is clearly seen from the right panel of Fig. 3, where the .ff -factor defined 
as K = daNLo/d&LO is shown. 

We have only considered QCD corrections to the gluon fusion subprocess in pp — > tiZ. 
Given the large gluon-gluon luminosity at the LHC and the fact that QCD corrections to 
gluon initiated processes are usually larger than those for quark-initiated, this seems to be a 
reasonable first step towards the complete NLO QCD computation. However, it is premature 
to discuss in detail how the NLO QCD corrections affect the analyses in Ref. 0,0]. We make 
one comment here. In Refs.0, H|, a(pp —>■ tiZ) was computed with the renormalization 
and factorization scales set equal to the top quark mass. As follows from Table [U for 
such renormalization and factorization scales the NLO QCD corrections for the gg — > tiZ 
subprocess are about 30 percent. Assuming that QCD corrections to quark intiated process 
are somewhat smaller, NLO QCD corrections of about 20% to pp — > ttZ for fi = m t can 
be expected. The NLO QCD corrections are positive and enhance the signal cross section, 
leading to a decrease in the luminosity required to probe the tiZ coupling. The NLO QCD 
computation also reduces the theoretical uncertainty of the signal cross section coming from 
uncalculated higher order corrections. We estimate that the residual scale uncertainty is 
±5% by varying the scale between /zq/4 and /i . 
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V. CONCLUSIONS 



In this paper we presented the NLO QCD corrections to gluon-initiated production of tiZ 
at the LHC This is a first step towards the full NLO QCD calculation of pp —>■ tiZ needed 
to study electroweak top quark couplings. For reasonable choices of the renormalization 
and factorization scales, the QCD corrections can be up to 75%. The remaining theoretical 
uncertainty in the cross section after inclusion of the NLO corrections is estimated to be 
±5%. The NLO QCD contributions are independent of the transverse momentum of the 
Z-boson. 

In the future, we plan to extend this computation to include the quark anti-quark annihi- 
lation and quark-gluon partonic scattering processes. It will then be possible to incorporate 
the improved knowledge of the signal processes into the analysis of Ref . j3, 0] • 
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